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Abstract 

Contribution to thermal conductivity from conduction electron scattering on crystal field mag¬ 
netic excitations is calculated and analyzed for normal rare-earth inter-metallic paramagnets. It 
is shown that in temperatures much lower than Debye temperature Tb its behavior essentially de¬ 
pends on the ground state of magnetic ion in crystal field and on the excitation energy in relation 
to Tb . Combined effect from the electron scattering on the crystal-field excitations, on acoustic 
phonons, and on nonmagnetic impurities is discussed in reference to CF splitting character and 
to the relative intensities of magnetic and non-magnetic scattering. Total thermal conductivity 
resulting from these three sources of scattering is calculated for REln^ (RE = Nd, Pr, Tm) and 
compared with experimental data. 

PACS numbers: 72.15.Eb, 72.10.Di, 71.70.Ch, 72.10.Fk 
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I. INTRODUCTION 


Crystalline field (CF) splitting of rare earth (RE) ions in metals influences conduction 
electron transport through their interactions with 4f electrons. In normal rare-earth systems 
(the ones with well localized f-level and a stable magnetic moment) Coulomb direct and ex¬ 
change interactions may be regarded as the most important. The last one gives magnetic 
contribution to transport coefficients which was studied theoretically already a long time 
ago. In early works concerning RE ions as impurities in metallic systems, it was consid¬ 
ered together with the contribution from isotropic Coulomb scattering.-^ The formula for 
electrical resistivity, derived in Ref. [2j] in frame of Boltzmann equation solution—and that 
for thermal conductivity derived in Ref. [lj in frame of the Kubo linear response theory— 
were then adapted for concentrated paramagnetic systems, respectively in Refs [s!^ 5]). The 
isotropic Coulomb scattering was omitted there, as being accounted in Bloch character of 


conduction electrons in a periodic lattice, and in Refs jj|5|] has been replaced by aspherical 
Coulomb scattering (electron-quadrupole scattering). The magnetic contributions to the 
electrical resistivity,- - and to the thermal conductivity^ have been expressed by a sum of 
components from conduction electron scattering on the 4f-electron magnetic excitations in 
crystalline held (inelastic - relating to transitions between the levels of different energies 
and elastic ones - between levels of the same energy). 

Calculation of the^ magnetic contribution to the electrical resistivity based on the for- 

erpretation of the experiment for 


6H8j. What seems more important, 


mula p mag (T), Ref. [3f] was successfully used in the int 
a number of RE paramagnets, see the examples in Refs 
some relations between the resistivity behavior and the character of CF-splitting have been 
established and experimentally confirmed^'-—. It was found that p mag (T) increases with tem¬ 
perature and in the range of temperatures corresponding to the energy of the first excited 
CF-level, the rate of increase depends on the excitation energy. In high temperatures the 
resistivity saturation value p spd , is independent of CF-splitting, while the zero temperature 
limit is governed by the 4f-electron ground state in the crystal held. Through these findings, 
the electrical resistivity behavior may serve as a kind of an identifier of the CF-splitting 
character. The situation is different in the case of the thermal conductivity studies. Despite 
it was found for some RE metallic compounds that the thermal conductivity experimental 
results may be approximated by its electronic part£2£ 5 the role of magnetic contribution in 
the electronic part behavior, as well as the behavior of magnetic contribution itself, is not 
clear yet. The formula of Ref. [.5] was applied to the low-temperature experiment interpre¬ 
tation for paramagnetic (or remaining paramagnetic down to a few Kelvins) compounds as 
TmSb, ErSb, PrPt 5 , PrCu 5 in Ref. p] for PrCu 5 in Ref. 11] or for PrAl 3 in Ref. 12]. All 
the considered compounds exhibit similar type of CF splitting - nonmagnetic ground state 
and comparable values of the first excited level energy. Their thermal conductivity behaves 
also in like manner so it is understandable that there had not been any discussion about the 
influence of CF-splitting character on the thermal conductivity temperature dependence. 
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This topic was also not approached in Ref. 13], where calculations of the magnetic contribu¬ 
tion to the thermal conductivity have been performed merely for two-level CF ferromagnetic 
system. However, the problem seems to be the current and relevant in view of the thermal 
conductivity experiment for REAl^ or /I'/Uri;;—. among which T 111 AI 2 , Prhi 3 , Ndhi 3 , or 
Tmhi 3 are paramagnetic or remain paramagnetic down to very low temperature. Their 
magnetic properties received attention in a number of works, including characteristics of 
their crystal field. The crystal field effects on the electrical resistivity of those compounds 
were also examined. Until now, however, there are no studies analyzing thoroughly their 
thermal conductivity behavior and the crystal held influence on it. 

We analyze this problem in the present work for RE inter-metallic paramagnets (RE— 
I-P). Their thermal conductivity we calculate as the effect of, additively treated, contribu¬ 
tions from the conduction electron scattering on CF-excitations, acoustic phonons and on 
the nonmagnetic impurities. Within Kohler variational approach,—^ we prove that, for the 
total thermal conductivity, and for their contributions from considered types of scattering, 


the simple formula can be used - the simpler one than that applied in Ref. 13]. For the 


electron-phonon contribution, we use the form after Kohler— and Ziman 1 -. For the impu¬ 
rity contribution, the variational formula—as we show —takes the standard form=^, derived 
within the relaxation time solution of the transport equation. The existence of relaxation 
time depending only on the electron energy we also justify for the electron scattering on 
the CF excitations. It ensures the Wiedemann-Franz law fulfilment in this case, analo¬ 
gously to the case of the electron-impurity scattering. In consequence, the formula for the 
magnetic contribution to the thermal conductivity, derived by us within the variational ap¬ 
proach proves as simple as that for the electrical resistivity derived within the relaxation 
time solution of transport equation in- 2 . It has simplified our analysis of the magnetic con¬ 
ductivity and allowed us to derive some rules governing the conductivity behavior in low 
temperatures as depending on the CF-splitting character. Combining this with the low tem¬ 
perature behavior of the electron-phonon and the electron-impurity contribution, we reach 
some general conclusions for qualitative behavior of the total thermal conductivity in this 
range of temperatures. Then, we verify these conclusions by comparison calculations with 
experiment for REhi :i (RE = Pr, Nd, Tm)—, and we discuss applicability of our model. Fi¬ 
nally, we refer critically to some aspect of calculation presented in Ref. 13(|, made within the 
similar physical model for the electron scattering and variational approach to the thermal 
conductivity computation. 


II. THE MODEL AND METHOD FOR CALCULATION 
OF THERMAL CONDUCTIVITY 

We consider below the contributions K a (T) to the thermal conductivity n(T) and p a (T) 
to the electrical resistivity p(T) from conduction electron scattering on acoustic phonons 
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(a = ph), on nonmagnetic impurities (a = imp) and on crystal field excitations (a = mag). 
Our considerations are based on the simple physical model, described in Appendices A-C. 

For the conductivity and the resistivity calculation we use Kohler’s formulas derived 
within variational method of Boltzmann equation solution^* 17 . We show (Appendix A) 
that for the applied physical model the simplest form of these formulas, (1A8D and (1A9D ) can 
be used. 

Applying in flA8l) - flA9D denotations P$ 2 = 3 P 22 j(Po{k B T) 2 7 r 2 ), P{\ = P^/Pq, and taking 
into account the form of scattering matrix elements P£{t), i = 1, 2, a = ph, mag (IBID . (1C3|) 
depending on the reduced temperature t — T/Td , we can represent the formulas in the form 


t 


k — K, 


Kn = 


Lq Td 

ry ? 

Po 


p a = PoPn(t) , 


Po 


TDOL 
r 0 




2 ’ 


( 1 ) 


where L 0 = n 2 k 2 / 3e 2 is the Lorentz number. 

The scattering matrix elements P°j, i,j = 1,2 for a = imp, mag we derive in Appendix B, 
and for a = ph we use the form (IBID . 

From the Matthiessen rule C(k, k') = )T« (7“(k, k') for the scattering probability (IA3D 
the corresponding rule follows for the scattering matrix elements and, consequently, for 
the total electrical resistivity p = J2aP a an d f° r fh e total thermal resistivity W(t) = 1 /n(t). 
For the reduced thermal resistivity W(t) = W(t)/Wo it takes the form 


W(t) = WP h (t) + W m ^(t) + PF imp , 

W^tj = W a (t)/W 0 , 

W a (t) = 1 /K a (t) , (2) 


where Wo = 1 /kq 1 , according to ([T]). 

Next, after substituting P 22 v (1B3D and P 22 ag (^) (EH1) in (JTj) we get from (TTT) — (O) 


W ph 



c mag _ 


ppmag _ c mag 



P r e _ p spd 

pS" ~ 


jmp 


w im p = 


c imp = 


ph" 

Po 


( 3 ) 


where P^ft) is described by (IBID and p res = Po mp Tii ip denotes the temperature-independent 
impurity part of electrical resistivity (residual resistivity) and p spd is the high temperature 
saturation value of p mag (t). 

The common way to find coefficients c mag and c imp (J3]) for particular RE-I-P relies on the 
assumption that the phonon contribution p ph (t) to the resistivity of this compound, p(t), 
and to the resistivity of its nonmagnetic counterpart, p NM (t), are the same. With the use of 
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Table I: Parameters characterizing relative impurity and magnetic contributions to the thermal 
resistivity. 


Compound 

J 

n ~ 1 ® ln/W 

£n=£ 1 

p res [pP cm] 

Pg h [pR cm] 

p spd [pR cm] 

Ref. 

c imp 

c mag 

Prlri3 

4 

0 

0.31 

7.5 

3.5 

20 

0.041 

0.022 

Ndln 3 

9/2 

1.2 

0.26 

8 

3 

j) 

0.03 

0.015 

Tmln 3 

6 

0 

1.8 

7.5 

3.1 

20.21 

0.24 

0.01 


approximations p ph (t 0 ) — 0 (to <C 1 ), p ph (too) = Po* too, (too 1 ) based on, correspondingly, 
the low and the high temperature asymptotics of p ph (t) (PQ), one can find p ph from 


Pnm (too ) = Pft oo + p 
Pnm(^-o) Pnm j 


res 

NM 


( 4 ) 


where denotes residual resistivity p imp (jT]) for the nonmagnetic compound. For magnetic 
compound, we use the low and the high temperature asymptotics of p mag (t) ( 0 Q) 


P mag (*o) = (pTVh) £ «i» , 

71=1 
£n=e l 

p^(too) = P spA = Po ag J(J + !)• 


( 5 ) 


Then, the residual resistivity for the magnetic compound p res and p ™ 8 can be found from 
the equations 


p(too') 

p(to) 


— p P O h too + p reS + p™ ag J( J + 1) , 


= P 


+ (pT g /h) 


'y ] n • 


( 6 ) 


Results of our calculation pQ h , p spd , p res and, next, c ma S 5 c im P with the use of dU)-© 
for RE I 113 (RE =Pr, Nd, Tm) are included in Table I. All the calculations were done 
with Tb = 170 K - following from specific heat measurements for Lalii 3 18 . To calculate 
J2 ™=i ain/h we have used CF energies and eigenfunctions from Ref. [lij corresponding to 

£n=£ 1 

x, W and J included in Table III. For all compounds we assumed t 0 = 5K/170K, and 
too = 300K/170K in correspondence to T D = 170 K and to the low temperature p(5K) and 
the high temperature p(300 K) resistivity values. These values we have read out from the 


electrical resistivity data in Refs [ 9 J20.21J and we have included them in Table II. For Prhi 3 
and for Ndln 3 we accepted as p^ the values p(5 K), respectively for Luln 3 — and for YI 113 -. 

In the case of Tmln 3 , we assumed that the compound has the same value of p^ h as 
L 11 I 113 , adopting for the last one the value following from Ref. 20]. We based on the remark 
that the high temperature slopes of the experimentally obtained electrical 


21 


(in Ref. 

resistivity graphs for the both compounds are the same. 
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Table II: Comparison of the electrical resistance value for REln ^, RE =Y, Pr, Nd, Tm, from 
various experimental works. (*) These data result from our extrapolation of the graph in Ref. |9] 
to T = 300 K 


Compound 

p(5K) [pRcm] 

p(300K)/p(5K) 

p(300K) [pRcm] 

Ref. 

Y Iris 

0.25 

28 

7 

_14 


0.33 

37.6 

13.2* 

.9 

Prln.3 

1.74 

3.5 

7 

14 


0.31 

48 

16 

20 

N dins 

1.46 

66.3 

17 

14 


1.5 

51 

16.5* 

9 

Tmlns 

1.8 

4.6 

8 

14 


1.8 

10 

18.1 

21 


In Table II we have also included the electrical resistivity values p (5 K), p(300 K) obtained 
from the experimental graphs for REln 3 (RE = Pr, Nd, Tm) and YIn 3 presented in Ref. fl3 ]. 
The puzzling differences between those values for P 1 T 113 , Tmln 3 YIn 3 , and the corresponding 




ones 


read out from Refs [9|j20|j21] we will discuss in Section IV. 


III. LOW TEMPERATURE CRYSTAL FIELD EFFECT 
ON THE THERMAL CONDUCTIVITY 

A. Crystal field effect on magnetic contribution to thermal resistivity 

We focus our attention on the range of temperatures 0.03 < t < 0.25, where the ther¬ 
mal conductivity for REI 113 (RE = Pr, Nd, Tm) 1 -, the RE- I-P representatives, shows pro¬ 
nounced qualitative differences. Examining the magnetic contribution to the thermal resis¬ 
tivity kP mag (J3J) in that range of temperature, we could confine our considerations to the 
scattering on the magnetic excitations between two lowest CF-levels. Henceforth, we take 
into account only the first component of both sums in P™ ag (1C3jl describing, correspond¬ 
ingly, the scattering on the inelastic and the elastic CF excitations. Further, we denote by 
an, 022 the sums of weights (aPj of all elastic excitations within the ground and the first 
excited level, and by of 2 - the sum of the weights of all inelastic excitations between these 
levels. Thus, we can express IU mag in a simple way, taking into account that two cases may 
be distinguished, depending on the ah value. 
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Table III: Characteristics of crystal field splitting for particular compounds, x and W denote the 
CF parameters, after given references. A is the energy of the excitation from the lowest level 
to the first excited, aR ( i,j = 1,2), defined in the text, characterize the weights of elastic and 
inelastic excitations corresponding to the lowest and the first excited level. For the characteristics 
we use CF energies and eigenfunctions presented in Ref. j3], correspondingly to given x, W and 
J. T d = 170 K. 


Compound 

CF characteristic 


X 

W [K] 

Ref. 

Ground and first 
excited CF state 

an 

022 

A [K] d = A/T D 

a 12 

Prhi3 

-0.66 

2.68 

22 

singlet-triplet 

0 

0.5 

101 0.6 

13.3 

Ndln 3 

0.31 

0.91 

23 

doublet-quartet 

6.7 

21.85 

6.9 0.038 

5.29 

Tmln3 

-0.665 

-0.57 

24 

singlet-triplet 

0 

0.5 

12.9 0.072 

27.97 


When an = 0 we deal with nonmagnetic CF ground state (NM) and denote IF mag = 
. Approximating Z(t) ~ 1 + exp (—d/t) (d = o?i 2 ) we get 


^NM g 

Wnm 

WSm 


c mag (W^M + W$ M ) , 
3 ai2 

t cosh 2 [d/2t] ’ 

3 Q>22 

t (exp ([d/t] + 1) 


(7) 


T^nm corresponds to the first sum in (IC3D . describing inelastic scattering, and WjjjM to 
the second one describing elastic scattering. 

When an ^ 0 (Zi > 1), the CF ground state is magnetic (M), and we denote IF mag = 
W™ e • Approximating Z(t) ~ Zi (1 + exp (—d/t)) we get similarly 


W™(t) 

W*(t) 


c^(ws +_K 1 ) • 

3 ai2 

lit exp [d/t] cosh 2 [d/2t\ 

3 an 

li t( 1 + exp [—d/t]) 


( 8 ) 


Values of ai] ( i,j = 1,2) can be calculated with the use of 4f-clectron eigenfunctions in 
the crystal field, presented in Ref. 19]. The ones calculated for considered here REln^ we 
include in Table III. As follows from the Table, the case of nonmagnetic CF ground state 
IRmag _ involves Prlii 3 and Tmlii 3 , while the case of magnetic CF ground state, 

W mag = W™ S i concerns Ndln 3 . 

We have analyzed (I71)-(lg|) for the considered range of temperature 0.03 < t < 0.25, and 
the range of d including values contained in Tabic III: 0.038 < d < 0.6. We present the 
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behavior of these functions (in said range of temperature), in dependence of d. in a schematic 
way in Table IV. Further conclusions follows from the relations 


Wnm q gi 2 W~m exp [d/t] an 

^NM a 22 2 <3,12 


( 9 ) 


fulfilled for the range 0.16 < d/t < 20 corresponding to the considered values of t and d. 

As it is seen from the Table III, in the case IF mag = IF/]/ 38 (ED, there is au 022 • It is 
the general property of the nonmagnetic CF ground state case, what may be verified with 
the use of 4f-electron eigenfunctions given in Ref.[19|. Hence IF^m 8 ma Y be approximated 
by its inelastic part — c mag VF ] ^ I , (ED, which is an increasing function for d > 0.3. For 

the energy d sufficiently small—with respect to lowest considered temperatures t > 0.03— 
the inelastic scattering becomes quasi-elastic and its contribution may be approximated by 
decreasing function IF^ 8 ~ 3c mag Oi2/F We have found this approximation being justified 
for d < 0.05. 

In the case IF mag = W™ e (JHJ) the weights of elastic ofT and inelastic ajT excitations 
may be comparable, as it is seen from the data for Ndlii 3 . For d sufficiently small, both 
parts Ibj^ and IF^ could give a comparable contribution to IF/!, 138 , which in this case may 
be approximated as W™ K = 3c mag (ai2 + 022 / 2 )/(lit). The component proportional to ajT 
describes the quasi-elastic scattering, whereas the one proportional to 022 describes the 
elastic scattering. 

For the magnetic ground state and large d (d > 0.3) the elastic scattering predominates, 
as it is seen from ( 03 ), and we get from flSJ) W™ e = 3c mag aH/(/i t). 

The conclusion is that for large d the predominant scattering is inelastic in the case 
of nonmagnetic ground state (NM) and is elastic in the magnetic one (M). For small d 
(d < 0.05) the predominant scattering is quasi-elastic in the NM case or elastic and quasi¬ 
elastic in the M case. Prlii 3 represents the first case, while Tmhi 3 and Ndlii 3 the third and 
the fourth cases, correspondingly. (For Tmlii 3 it concerns t > 0.05, due to d value, see 


Table III.) 

Noting that the temperature dependence of PF mag for the three last cases (elastic, quasi¬ 
elastic, elastic plus quasi-elastic) differ only by coefficients, we identify all them and describe, 
by contract, as elastic scattering. We denote for these cases IF mag = IF™ 38 . 

Using in H /7 ^ n M (J7|) the approximation cosh(d/2t) 2 ~ exp (d/t) justified for large d, one 
can note that IF,™ 38 is about exp (d/t) times smaller than IF™ 38 ~ 1/t, accurate to the 
corresponding excitations weights ratio. 

Identifying d > 0.3 as a large energy, and d < 0.05 as a small one, we can summarize: 

(i) If the electron scattering on CF excitations is predominantly inelastic (nonmagnetic 

RE ion ground state and large excitation energy d) the magnetic contribution to the thermal 
resistivity IF mag = IFF 38 has the form (JTJ) , increasingly depending on temperature. 

(ii) If the CF-scattering is predominantly elastic (magnetic ground state and large ex¬ 
citation energy d) or quasi-elastic (nonmagnetic ground state and small d) or elastic plus 























Table IV : Type of monotonicity for inelastic and elastic parts of magnetic contribution to the ther¬ 
mal resistivity IV ma §(f) © approximated for 0.03 < t < 0.25 for magnetic {Wy\) and nonmagnetic 
(IVnm) ground CF state. Symbols Z\ \ denote respectively the function increasing or decreasing 
in the whole interval. Description max denotes maximum of the function within the interval. 



Wnm 




d > 0.3 

Z 1 

/' 


\ 

0.05 < d < 0.3 

max 

max 

max 

\ 

d < 0.05 

\ 

\ 

\ 

\ 


quasi-elastic (magnetic ground state and small d) the magnetic contribution to the thermal 
resistivity W mag = is decreasing function. For all these cases it can be approximated 

W™ g ~ a/t where a = 3c mag on/Z in the hrst case, a = 3c mag ai2 in the second case, and 
a = 3c mag (oi2 + 0^22 /2)/Zi in the third one. 

(hi) W™*~eM-d/t]W^. 

In Fig. 1-Fig. 3 we show VF mag Q calculated for all crystal held levels corresponding 
to CF splitting for Prlii 3 , Ndlii 3 , and Tmhi 3 , (according to x, W given in Table III). As 
c mag we use the values included in Table I. The graphs IF mag confirm (i)-(iii), however for 
Tmln 3 (Fig. 3) it concerns t > 0.05. Only for this range of temperature the CF-scattering 
for this compound is quasi-elastic and Tmln 3 may be included into the same group (ii), 
to which also Ndhi 3 belongs. For both compounds W mas = W™ g is approximated by 
the same formula IR mag ~ W e “ ldg = a/t , where a = 3c mag ol2 corresponds to Tmln 3 , whereas 
a = 3c mag (012 + 022 / 2 )//! corresponds to Ndln 3 . Calculating a for both compounds, with the 
use of data from Table I and Table III, we get the greater value for Tmln 3 . The same relation 
holds for the values of IV mag for these compounds, what is seen comparing corresponding 
graphs in Fig. 2 and Fig. 3. We can recognize this as an evidence of the correctness of our 
low-temperature approximation IF e “ iag ~ a/t in case (ii), at least with respect to Tmln 3 and 
Ndln 3 . 


B. Crystal field effect on total thermal conductivity 


The low temperature behavior of the thermal conductivity for considered here paramag¬ 
netic REln 3 (RE = Pr, Tm, Nd) 14 , is similar to the behavior observed for their nonmagnetic 
counterparts YIn 3 and LUI 113 . The conductivity of Prhi 3 and Ndhi 3 shows maximum like 
that of Yln^ while the one of Trains increases like that of LUI 113 — . 

For nonmagnetic metals, both types of behavior have long been known and were explained 
in Ref. 26] within the employed here model. Thus, we can repeat the arguments presented 
there with the use of K,(t) = l/IF(t), ([2]), after omitting magnetic scattering. Approximating 
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W ph (t) (ED by its low-temperature form W ph (t —t 0) we get 

-i -i 


K(t) = iyp h + iy im p 


fyp* 1 ~ 12 J 5 (oo)n s 


t 2 


71 


( 10 ) 


where J$(oo) ~ 125. (Js(t) and n s are defined in (IB1I) ). Behavior of «(£) (TTOT) is the effect 
of a competition between inelastic scattering on phonons - described by increasing function 
W ph (TTUT) and elastic scattering on impurities - described by decreasing function W im p (ED- 
It can be easily verified that n(t) (1IU1) exhibits maximum at 


7T 2 c im P \ 1/3 

3 000 n s J 


( 11 ) 


Requirement that the maximum occurs inside the considered interval of temperatures, 
0.03 < t max < 0.25, gives the upper and the lower limit for parameter c imp /n s . For a fixed 
n s , there is a threshold value for c imp = p res /Po h , beyond which the maximum disappears in 
this range of temperatures, and k(£) becomes an increasing function. For fixed p ph , we can 
consider the threshold value for p res . 

A direct conclusion from 03 is that hi(t) adopts greater values for smaller c imp or 
smaller n s . Applying numerical analysis, we have found that the maximum becomes more 
pronounced for smaller c imp , but flattens for smaller n s . 

Consider now the low temperature behavior of the thermal conductivity for magnetic 
compounds - described by the formula n(t) = 1 /W[t) (ED • Recalling (i)—(ii) from the previ¬ 
ous subsection, we can note that behaves similarly to W ph (t ) (fTOl) and behaves 

in the same way as W imp . It is also seen, in Fig. 1-Fig. 3 for PF ma s = IF™ e f j s corresponding 
to Prlii 3 and for W ma,g = IFg™ 518 corresponding to Ndln 3 , Tmln 3 , where graph of W ph (t) 
represents the non-approximated form of the phonon contribution (ED- Combining above 
comments with the theory for nonmagnetic metals^, we get its extension for magnetic met¬ 
als (RE- I-P). 

(I) In systems with a predominance of inelastic CF scattering (see (i) in Section III A), the 
thermal conductivity exhibits the low temperature maximum when the impurity contribution 
to the scattering IR im P is sufficiently small with respect to the contribution from the sum 
of and • 

(II) In systems with a predominance of elastic or quasi-elastic (or elastic plus quasi¬ 

elastic) CF scattering (see (ii) in Section III A), the thermal conductivity exhibits maximum 
when the sum of contributions from IR im P and is sufficiently small with respect to 

W ph . 

Realizing that acts like reducing while acts like enlarging W imp in nonmag¬ 

netic metals, it is easy to see that £ max - the temperature of the maximum and its value 
depend on c imp in the similar way like for nonmagnetic metals. Therefore, in each case (i) 
and (ii), for a fixed n s and c mag there is a sufficiently small value of c imp required for the 
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Figure 1: Reduced thermal resistivity W(€)/Wq = W(t) ((2j) and its contributions W a (t) /Wq = 
W a (t) fl3J) for Prln 3 . n s = 1.3, c mag = 0.022, c imp = 0.15. 


conductivity maximum in the range of temperatures 0.03 < t < 0.25. Similarly, there 
must be a threshold value for c imp = p res /Po h (a threshold value for p res when Pq 11 is fixed), 
beyond which the maximum disappears in the considered range of temperatures, and the 
conductivity becomes an increasing function. 

In contradistinction to the nonmagnetic metals case, we cannot establish its value ana¬ 
lytically. However, for RE- I-P of the same phonon contribution (fixed n s and pg h ) some 
conclusions follow from (I)—(II) combined with (i)-(iii): 

(Rl) In case (I), the threshold for c imp is the larger the greater is (the greater is 

c mag a 12 and smaller is d). In case (II) the threshold is the larger the smaller is (the 

smaller is a). 

(R2) The threshold in case (II) should be smaller than that in case (I), independently of 
values of c mag corresponding to IF™' 18 and IF™)) 8 in these cases. 
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t=T/T D 


Figure 2: Reduced thermal resistivity W(t)/Wo = W(t) ((2j) and its contributions W a (t)/Wo = 
W a (t) (O for Ndln 3 . n s = 1.3, c mag = 0.015, c imp = 0.08. 


Applying (Rl) to Ndln 3 and Tmln 3 we can conclude that the threshold for the first 
compound should be greater, because W e “ ag is greater for Tmln 3 . It follows from the 
greater value of a for the last compound, which we have noted in the previous subsection, 
and what is seen from comparison of the graphs bF mag (f) in Fig. 2 and Fig. 3. Applying 
(R2) to RE In 3 (RE =Pr, Nd, Tm) we may expect the threshold for Prln 3 being greater 
than for the other compounds. Our numerical findings confirm these relations, although the 
difference between thresholds found for Ndln 3 and Tmln 3 is slight. The obvious conclusion 
from relations between thresholds is that the conductivity of Prln 3 may exhibit the more 
pronounced maximum for the value of Ci mp (p res ) greater than a value corresponding to the 
conductivity maximum for Ndln 3 . On the other hand, the value c imp (p res ) for which the 
conductivity of Tmln 3 increases may be sufficiently small for the conductivity of Prln 3 to 
exhibit a maximum. 
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Figure 3: Reduced thermal resistivity W(t)/Wo = VF(t) <[2j) and its contributions IF a (i)/H o = 
W°(t) © for Tmln 3 . n s = 1.3, c mag = 0.01, c imp = 0.15. 


The above conclusion is confirmed by the conductivity experiment for REln^ ( RE = 
Pr, Nd, Tin) 1 -, when one takes into account the residual resistivity values p res (NdIn 3 ) = 
0.26/rflcm, p res (PrIn 3 ) = 1.74/iO cm, p res (TmIn 3 ) = 1.8/iO cm (obtained from (© with 
p(t 0 ) = p(5K) included in Table II, p“ ag from Table I and «=i a\ n /li corresponding to 

e n =e l 

CF splitting data from Table III. 

We also illustrate this with the calculation of W(t) (1211 . what is seen in Fig. 1-Fig. 3 when 
one takes into account the equivalence between the conductivity maximum and the resistivity 
minimum, and between the conductivity increase and the resistivity decrease regions. 

As the additional difference between the conductivity for Prln 3 and for both Ndln 3 and 
Tmln 3 we can note its much greater values than for the other compounds. It is seen from the 
experiment data presented in Fig. 4, and also from results of our calculations of W(t ) (J2J) in 
Fig. 1-Fig. 3. This is a simple consequence of relation (iii) in Section III A: W™\ g !Fe} ias . 
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Figure 4: Thermal conductivity for Prln 3 , Ndln 3 , and Tmln 3 aftei— . 


We can generalize these characteristics of the crystal held influence on the conductivity 
for RE ln 3 (RE =Pr,Nd, Tm) for arbitrary group of RE- I-P of the same W ph . The gen¬ 
eral conclusion appropriate for such compounds is that the conductivity of compounds of 
type (i) is more likely to exhibit a maximum and larger values than the conductivity of the 
compounds of type (ii). 


IV. COMPARISON WITH EXPERIMENT 

With the use of the theory presented in the previous section and the results of calculation 
of W in Fig. 1-Fig. 3 we explain the qualitative differences between the thermal conductivity 
experiment for RE I 113 (RE =Pr, Nd, Tm)—, likewise the relations between the conductiv¬ 
ity values. These results, however, cannot be used for a thorough comparison with the 
experiment—. 
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t-T/T D 

Figure 5: Thermal conductivity n(t) = 1 / (Wq W{t)) (J2|) for CF splitting corresponding to Prln 3 , 
Ndln 3 and Tmln 3 described in Table III, c mag , c imp from Table I, and n s = 1.3. 

K ph = 52.1 [W/mK] corresponds to pg h = 8 [//Hem]. 


.P h 


r> ma g 


For this purpose the calculations should be done with c imp = p res /po 
P spd /P 0 h J(J + 1) following from the electrical resistivity experiment for REln 3 , RE = 
Pr,Nd, Tm, as well as for their nonmagnetic counterpart YIn 3 , which was performed with 
the use of the same samples as were used for the conductivity results in Ref. [bij], and was 
presented therein. Following the approach described at the end of Section II, we have got 
p(t 0 ) and p(too) reading out the low temperature p(5K) and the high temperature p(300K) 
resistivity values from the graphs^. As p^|m we have accepted the value p(5K) for YIn 3 . 
(All these values are included in Table II). However, substituting p^M and p(to), p(foo) for 
Prhi 3 and Tmh. 13 , to (J3D—(JHD we have ended up in contradictions. Therefore the method 
which works_for the resistivity data obtained from Refs 20 2l|, fails for the data obtained 
from Ref. [l^]- If one regards that in the last paper p(300K) 


or Prln 3 , Tmln^ , and YIn : 


3 , 


3 . 2 O. 2 I] (see Table III), the 


is about two times smaller than the respective value in Refs 
failure of the method is not such surprising as the results of Ref. [14j themselves. 
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Comparing the residual-resistance ratios (RRR) p(300 K) / p (5 K) and the residual- 
resistance values p res ~ p(5 K) included in Table II, one can notice that, for Prln 3 , the 
results presented in Ref. 14] are_obtained with the use of a sample of much lower purity 
than the one presented in Ref. 20]. It may be an explanation of the great differences be¬ 
tween the high temperature values of the resistivity for this compound obtained in Ref. 14] 


and in Ref. 20]. As it was argued in Ref. 17], a large contribution of residual resistivity to 


the total resistivity essentially influences the character of electron-phonon scattering and 
makes this scattering dependent on the electron-impurity scattering. For such a case, our 
way of calculation based on the Matthiessen rule, including equations (J!])-©), may be inap¬ 
propriate. Such explanation, however, does not fit to the results for Tmln 3 and YIn 3 . The 




20 


RRR values for these compounds^ are about half of the values following from Refs [9 
whereas the values p res ~ p(5K) 14 are comparable with the corresponding ones in the first 
two papers. We must recognize the differences between the results for Tmln 3 and YIn 3 in 
Ref. 14j and in Refs @1 2o| as unexplained. There is a different case with Ndln 3 , for which— 


14 


agree 


in contrast to the other compounds—all the electrical resistivity values in Ref. 
very well with those in Ref. [9]. It justified, in our opinion, accepting the values c imp , c mag in 
Table I as corresponding also to the sample of Ndln 3 used in Ref. 14]. Therefore, among 
the results of our calculation of n(t) = l/(Wo W[t)) (J2]) presented in Fig. 5, the only reliable 
comparison with the experiment— may be performed for Ndln 3 . Comparing graph of n(t) 
for this compound in Fig. 5 with the experimental behavior of thermal conductivity in Fig. 4, 
and taking into account the temperature scaling t = T/T b with T f> = 170 K one can note 
that: 

- K,(t) in Fig. 5 has the maximum at temperature about twice smaller than the experi¬ 
mental one Lav — 0.17, and for the temperature of the minimum t min ~ 0.37 the difference 
is greater. 

- values of n(t) are smaller than the corresponding experimental ones, and this dis¬ 
agreement, being about 20% at the temperature of n(t) maximum, reaches 30% for other 
temperatures. 

It should be noted that we did not obtain better agreement for other values of n s . We 
have found, similarly as in the case of non-magnetic compounds discussed in Section III A, 
that increase of n s moves the maximum towards lower temperatures, and that decrease of 
n s cause its flattening. 

There are several possible causes of discrepancies described above. As the crucial one, 
we consider deficiencies of the standard model for the electron-phonon scattering and the 
electron scattering on nonmagnetic impurities, which we used in our calculation. 

In the simplest way it may be demonstrated by comparing n(t) (HUD with the experiment 
for YIn 3 —. The thermal conductivity for this compound shows the maximum at t max = 
0.17 assuming Tb = 170 K, as we have assumed in the calculations of previous section. 
Substituting this in (TTIT) together with c imp ~ 0.06 - following from the resistivity data in 
Ref. [l4|, or together with c imp ~ 0.03 - following from the data in Ref. [jj, one obtains 
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the nonmetallic values of n s . On the other hand, for metallic values (n s > 0.63) and 
for 0.03 < c imp < 0.06 one gets 7R > 400 K which seems to be unrealistic. The general 
arguments for the theory oversimplification can be found in Refs [131 27, Q. They refer to the 
electronic transport in non-magnetic metals, but they remain relevant also for the magnetic 
ones considered in the present paper. The unclear sense of parameter n s with respect to 


these metals, which we have mentioned in Ref. 29] is the one example. Another one is the 
application of the Matthiessen rule. We have recognized it as inappropriate for calculation of 
thermal conductivity in case of low purity sample of Prlii 3 in Ref |l4|], but its application is an 
approximation in every case, and leads to a discrepancy with the experiment. We may expect 
this discrepancy increase due to the scattering on magnetic CF excitations. As follows from 
the thermal conductivity calculation for the CF-split RE impurities in a metal-, there are 
essential differences in the thermal conductivity values depending whether the Matthiessen 
rule is applied or not. 

As the additional cause of poor agreement with the experiment in the case of Ndln 3 
we may consider the omitting—in our calculation—the aspherical Coulomb interaction as 
a source of electron scattering. The important role of that scattering for the electrical 
resistivity for Ndln 3 has been proven in Ref. o|. 

In the case of Prln 3 and Tmln 3 , although the results of our calculation presented in Fig. 5 
have no direct relevance to the experiment—, they agree qualitatively with that experiment. 
It concerns both the conductivity behavior for particular compound as well as relations 
between the conductivity values for all REln 3 (RE = Pr, Nd, Tm). We can conclude that 
the CF-scattering effects in these compounds are strong enough to make their conductivities 
less sensitive to the effects of the non-additivity of the scattering intensities. 


V. SUMMARY AND CONCLUSIONS 

We have analyzed how the conduction electron scattering on magnetic crystal-held (CF) 
excitations influences the low temperature dependence of the thermal conductivity in rare- 
earth intermetallic paramagnets (JLF-I-P). We have considered the thermal conductivity as 
resulting from the conduction electrons independently scattered by acoustic phonons, by the 
non-magnetic impurities and by the CF magnetic excitations. Hence we have applied the 
formula hi = l/{W ph + W imp + W ma9 ), where W° the a = ph, imp, mag are contributions 
to the thermal resistivity W = 1/n. They are proportional to the material parameters p ph , 
p res , p mag , which can be found from the electrical resistivity experiment in a way based on 
the Matthiessen rule, what we have described in Section I. 

Calculating the magnetic contribution and examining its behavior below 0.25 To , we have 
found that it increases like W ph when the magnetic scattering is dominated by the inelastic 
CF excitations and decreases like W imp when the elastic and/or quasi-elastic scattering 
prevails. In the first case (i) we denote W ma9 = IF™) 5 , in the second one (ii) W ma9 = W ™ a9 , 
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considering also the quasi-elastic scattering as the elastic one. We have found that <C 

W ™ a9 and established dependence of W ™ 9 and W™ a ° on some parameters characterizing 
CF splitting. The case (i) concerns compounds like Prln 3 , for which the f-electron ground 
state is nonmagnetic and d, the energy of the first excited CF state, is sufficiently large in 
comparison with 0.25 Tb • The case (ii) concerns compounds where d is sufficiently small, 
like Ndln 3 (with magnetic ground state) or Tmln 3 (with nonmagnetic ground state). We 
illustrated this in Fig. 1-Fig. 3 calculating W ma9 = W ma9 /W 0 (JSiH©- 

From these hirelings about W ma9 we could conclude that values and behavior of the total 
thermal conductivity in RE- I-P, like those in nonmagnetic metals, results from a competi¬ 
tion between contributions from inelastic and elastic scattering. The inelastic scattering is 
described by W ph and in case (i) additionally by W ™“ 9 , while the elastic scattering is de¬ 
scribed by W imp and in the case (ii), additionally, by W ™ a9 . When the elastic scattering is 
sufficiently small with respect to the inelastic scattering, the conductivity exhibits low tem¬ 
perature maximum, otherwise it increases. For a particular RE- I-P, like for nonmagnetic 
metals, the temperature of the maximum increasingly depends on c imp , or, equivalently, on 
the value of p res - both determining the impurity contribution. However, it depends not 
only on W ph , as in the case of nonmagnetic metals, but also on W ma9 . The same concerns 
the threshold for c imp (p res ) , beyond which the maximum disappears. The dependencies of 
the thresholds on parameters characterizing CF-splitting we described in (Rl) Section 111 
for both cases (i) and (ii). We have also established (in (R2) therein) that the threshold 
value in case (i) is greater than in case (ii) and that the same relation concerns the con¬ 
ductivity values. These relations, which apply to RE- I-P of the same phonon contribution, 
lead to the conclusion that the conductivity maximum and its large values require much 
purer samples in the case of predominately elastic and/or quasi-elastic CF-scattering than 
in the case of the inelastic scattering predominance. It can be seen analyzing the conduc¬ 
tivity experiment for REI 113 ( RE = Pr, Nd, Tin) 1 -, and corresponding values of p res , as we 
have discussed in Section III. We have also illustrated it by calculating the reduced thermal 
resistivity W(t) for these compounds, presented in Fig. 1-Fig. 3. 

In these calculations we could not use the values c imp = p res /p ph corresponding to the ex¬ 
perimental data in Ref. [14], because hireling a consistent set of material parameters p res , p ph , 
p mag , based on the electrical resistivity results in that paper for RE In 3 (RE =Pr,Nd, Tm) 
and for their nonmagnetic counterpart YIn 3 , proved impossible. As we note for Prln 3 , 
Tmln 3 , and YIn 3 the resistivity experimental values in Ref. 14] differ greatly (in the whole 
investigated range of temperatures) from the corresponding ones in Refs fol 2ol 21]. In case 
of Prln 3 , basing on the residual-resistance ratio (RRR) and p res for the sample used in 
the experiment^, we could ascribe this difference to the Matthiessen rule breaking caused 
by low purity of the sample. For Tmln 3 and YIn 3 we had to recognize the differences as 
unexplained. 

The consistent set of material parameters and, consequently, the coefficients c mag and 
c imp for RE In :i (RE =Pr, Nd, Tm) could be found in the electrical resistivity experiment 
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of Refs {9 


20 


21]. These coefficients (included in Table 1) we used in our calculation of 
n(t) = 1/(WoW{t)) (ED-© and presented results in Fig. 5. The only reliable comparison 
of these results with the experiment— we could do for Ndlii 3 , for which the values of the 


electrical resistivity in Ref. 14] agree with those in Ref. |9]. From discussion of discrepancies 
between the results of our calculation in this case and the experiment^ one can conclude 
that our model of calculation is too simple to provide a good quantitative agreement. The 
obvious way of improving the model is to go beyond Mathiessen rule in calculation, what 
seems to be a nontrivial problem to solve. The other improvement may be to take into 
account the aspherical Coulomb scattering as the additional source of scattering, but it 
seems to be of specific (for Ndlri 3 ), not of a general importance for RE- I-P. Nevertheless, 
the results of calculation presented in Fig. 5, also for Prlii 3 and Tmln 3 —although their 
quantitative agreement with the experiment^ is worse than that for Ndln 3 —incline us to 
believe that the model may be useful for qualitative interpretation of thermal conductivity of 
any RE- I-P. It proves to be particularly useful for analyzing the CF effects on the differences 
in the conductivity behavior for a group of iso-structural RE- I-P, as we have shown in our 
discussion performed for REln 3 . 

It should be stressed that all these capabilities of the model are based on the simplicity 
of the variational formula for the conductivity (IA8j) . much simpler than that was used in 
the thermal conductivity calculation for ferromagnetic RE intermetallics in Ref. 13]. As 
we have shown in Appendix B, the term — T S 2 Lee, of which the latter formula differs 


from 


is negligible for the physical model we have considered. It can 


:>e similarly 


demonstrated that this term is also negligible for the model considered in Ref. 13]. In that 
paper it was calculated with the use of incorrect form of the electron-phonon contribution 

31| . It caused, as it seems, its 


30 


to the thermoelectrical power, S ph obtained in Refs 
overestimation. We have proved incorrectness of the formula for S ph — in Ref. 32]. 


Appendix A 


To get the formula for thermal k and electrical a conductivities it is usual to start from 
linear relations between the electrical and the thermal currents on the one hand, and the 
electrical field and the temperature gradient on the other hand, with the coefficients Let , 
Lte , Lee , Ltt ,( the generalized transport coefficients), see Ref. [13], Chapt. VII. With the 
use of expression S = — L ET /L EE for the thermoelectric power one can write 


k — —T S 2 Lee — Ltt , 

o' = Lee ■ (At) 


The transport coefficients can be found with the use of the Kohler’s variational method 
for the Boltzmann equation^ - —. Its solution, the non-equilibrium distribution function, 
is represented there as the linear combination of some basis 0j(k), i = 1 and the 
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coefficients of this combination are found from the so called variational principle. After 
decomposition of the distribution function with respect to the basis one gets the generalized 
transport coefficients expressed by trial currents ,/,, [/* 


Ji = 


d k 


-df? 


Uj = —e dk 


de k 

'zdf? 

de. 


k J Mk)(v-u) , 

‘ A.(fc)(*> ■«)(£*-C (T)), 


(A2) 


and the elements of the scattering matrix P,j, i,j = 1 ,,n: 

V 


p ij(T) = -pf /dk/dVC(*,fc')/»(l-/“) Mtf (fc,V) 


(A3) 


Uij{k, k’) = [c fa(k ) - 0j(fc')][^i( fc ) - ')]> 

where f k is the electronic equilibrium distribution function, u denotes the unit vector in 
the external field direction, v is the electron velocity, e k = ( hk) 2 /2m - its energy, ( (T) - 
the chemical potential. C(k, k') denotes the transition probability per unit time for the free 
electron scattered from the state k to the state k'. 

Considering electron scattering on phonons, nonmagnetic impurities and CF-excitations 
we will use denotation C a (k,k') (a = ph, imp, mag), and P“(T) for the corresponding 
scattering matrix element (lA3jh In this way, for each a we can specify generalized transport 
coefficients L^ u (/i, v — E,T) and similarly n a , cr a , S a - the a contributions to the thermal 
conductivity, the electrical conductivity, and the thermoelectric power, respectively. The 
expressions for the generalized transport coefficients in the variational approximation of the 
n-th order were derived in Ref. 15] with the use of the base functions 

0i(fc) = (k-u) [e k - CCOP 1 , % = 1,..., n. (A4) 

In the second order, the transport coefficients Tg E , T , for every a have the form 


J EE 


7 TT 


= JL ( 1 , Jj Pi _ 2 J_i Er 

j i P22 Ji P22. 


Pi 


11 


(i + El 

T P 22 V Ui P u 


JhPu 

U 2 Pn 


(A5) 


where—for simplicity—we have omitted symbol a in L“ v (fi, v — E,T) and in P“(T). 

To estimate relations between the components of the sum in (1A5D we use trial currents 


relations J 2 /J 1 = 7r 2 (k B T) 2 /e F , U\/U 2 = 3/(2£ F h following from their form derived with 
the basis (1A4D i — 1,2, see Eq. (9.12.11) in Ref. [17]. Taking into account additionally the 
relations between the scattering matrix elements, following from their forms (IBID . (IB3D . 
and (IC3D it is easy to verify that, for every a = imp, ph, mag, the second and the third 
components in parentheses in (1A5|) are negligible in comparison to the first one. Hence the 
following approximations are justified for every a 


Lee — 


Jl 


P^ 


T'P'P ~ — — 


Pf 


11 


T P 


22 


(A6) 
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where J\ = e k F /{3n 2 h), U 2 = Ji(vr k B T) 2 /(3e),—, and is the Fermi radius. 

Now we compare the first component of k (I A1 [1 . with the second one (IA6[) considering 

relation T S' 2 L EE /|L TT |. 

For a = mag, imp, we use the form S a = —7r 2 k B (k B T)/(3e£ F ), which can be derived, 
with the scattering matrix elements (1B3D in the case a = imp and the elements (1C3D in 
the case a = mag, in the same way as it has been done in Ref. 32] for the conduction 
electron scattering on the magnetic levels of 4f-electrons in the molecular held. Substituting 
in (1A6j) the scattering matrix elements (1B3D and (1C3D we get for the both considered cases 
TL ee /|L tt | = 3e 2 /7r 2 kl and finally TS 2 L ee /\L tt \ = K 2 (k B T) 2 /3{£ F ) 2 < 1. 

In the case of a = ph, we use (S' 1311 ) 2 < 7r 4 k 2 (k B T) 2 /(£ F ) 2 , following from the form S ph in 
Ref. 17], (9.12.20). With the use of the scattering matrix elements (IB 1|) and (1A6|) we get 


S 2 TL 


EE 


< 


(k B Tf 


7 TT 


3 7T + 


9 n s 
t 2 


= R ph 


(A7) 


Because R ph < ( k B T) 2 (3n 2 + 9 n s )/£ 2 for t > 1 and R ph < (k B T E ) 2 (3n 2 + §n s )/£ 2 for 
t < 1 we obtain that TS 2 L ee /\Ltt\ 1 for every t. Summarizing, for every a and all 
experimentally accessible temperatures, one can neglect the first term in the r.h.s. of (1A1[) 
and approximate n a (T) ££ —L^ r . Obeying (IA6I) we get 


n a (T ) 


vr A Jl (lk B Tf 
9 e 2 T P? 2 ' 


(AS) 


Similarly, substituting L EE (1A6ji into (1A1I) we get the a contribution to the electrical 
resistivity p a = l/a a 



(A9) 


Appendix B 


To deal with acoustic phonons we assume Debye model and deformation potential 
approximation—. The scattering matrix elements P? h (1A3D derived with the use of the 
basis (1A4I) i = 1,2 can be found in Refs 15-17]. In terms of reduced temperature t = T/T E 


Pfi =AP ph t 5 J 5 (l/t), 


pph _ pph 

- r 12 — y r ll ) 


it =T ( KT ) 2 Pf? 


where 


dz 


_ 3 n s 

1 H- f - 

7r 2 t 2 


1 


sinh 2 (z/2) 


2vr 2 

= ^Jn(x) , 


(Bl) 


are Debye integrals, n s = k 2 F /q B 2 is defined by the Fermi k F and the Debye q D radii. £ s = 
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2 m v 2 is the energy of electron of the wave vector q s = 2mv s /h, where v s is the sound 
velocity averaged over directions in a crystal. Value of pP h = 2 C 2 Vom 2 g D 5 /(247r 3 ^, 4 Mv s ) is 
expressed by C — the energy of electron-phonon interaction, Vo - volume of primitive cell, 
m - mass of free electron, M - mass of an ion. 

Calculating scattering matrix elements P™ p with the base (1A4|) we assumed the standard 
form of the scattering probability for the conduction electron potential-scattering on the 
ionized impurities,— 


C imp (k',k) = C^ p g imp (q)S(£' - e) , 
imp 32 (tt) 3 ZV rP^P 

0 V he 2 A 4 ’ 

9 {q> (1 + (?/A) 2 ) 2 ' 


(B2) 


where A is the screening constant, n imp - the density of the impurities, Z - the effective 
valence of the impurity, and e - the dielectric constant. In the standard approximation™, 
A ~ 8 Cf/3. 

Performing the integrals in (IA3I) we use identity f°(ek)( 1— f°(e = k B T {—df(£k)/d£k) 


and then apply the Sommerfeld expansion (see e.g. Ref. 33]) confining to the first non¬ 
vanishing term in the approximation of the strong degeneration of the electron gas. 

We obtain 


where 


pi m p 


pimp 
9 22 


Po mp H^hp/X) 
( k B T ) 2 ^ 


7r 


(. k,r > 2 p;r 


Hv) 

phn p 


ln(l + p 2 ) 

4 Z 2 e 4 
37 t e 2 
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1 +p ’ 


m 

h 5 


n 


imp 


(S3) 


(B4) 


The form (1B2D of scattering probability ensures the existence of the relaxation time 
solution for the Boltzmann equation. 31 That is why, substituting (1B3D in (IA8[) and in (1A9|) 
one gets the standard results for the impurity contribution to the thermal conductivity n, imp , 
and for the impurity contribution to the electrical resistivity (the residual resistivity p res ) 
17 26 . These contributions obey the Wiedemann-Franz law which is the consequence of the 
relaxation time existence, the same for the charge and the heat transport. 
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Appendix C 


The scattering probability for the conduction electron scattered on the CF magnetic 
excitations in a cubic crystal has the form— 

2J+1 

£'-£) + E cn 

71 = 1 

Em>E n 

cr = Y (iex^-1)) 2 ^, (Cl) 

where j ex denotes the energy of s—f exchange interaction, and g is the Lande factor. The first 
component in the external brackets refers to the elastic scattering on the zero-frequency CF 
excitations, the second one to the inelastic scattering on the CF-excitations of the energies 
l^rnn — i(F/ m E n ) 7 ^ 0 (^rim A nm /k bT) and 

2J+1 

9o ~ 3 a nm p n , 

71=1 

En — Em 

gZ s (T) = 3aZiPn - Pm ), = |<n| J»| 2 , 


(C2) 


C' mag (k', k) = 


Cg ias 





where J denotes the total angular momentum of 4f electrons, |n), | m) are the 4f eigen-states 
in the crystalline field and E n , E m are the corresponding energies; l n is the degeneracy of 
the level of the energy E n , and p n = exp[— E n /k B T]/Z(T) is the probability of this level 
occupation (Z(T) = J2n ^ exp [—E n /k B T]). J z is the z component of the operator of total 
angular momentum. The matrix elements a zz m have the meaning of the excitations weights. 

We will calculate the scattering matrix elements P™ ag i,j = 1,2 (1A - > ll with the use of the 
base functions (1A4|) in the way described in Ref. 32], based on the property of C' mag (k',k) 
being even function with k and with k'. This property, similarly as in the case of electron- 
impurity scattering ensures the existence of the relaxation time solution of the transport 
equation—, the same for the charge and heat transport. We write the final form of F*™ 3,8 , 
i,j = 1, 2 as depending on reduced temperature t = T/T D 


Pn S (t) 


pmag 


p mag 
1 2‘2 


pmag 


pmag 


2J+1 

£3 


71=1 


1 - exp [~d nm /t\ 
sinh [d nm /t] 


£m >£n 


2J+1 

+ £ 3 CLmnPn, 

71=1 

£ n =£ m 


2 pmag 

-(F b T) 2 ^ 

O £t? 


y(fc B T) 2 Pr g 

2Vo3ex(9 ~ l) 2 m 2 K 
37r 3 h 5 


(C3) 


where e n = E n /k B T D , e m = E m /k B T D , d nm = A nm /k B T D , and p n = exp [-e n /t\/Z(t) with 
Z(t) = l n exp[-£ n /t]. 
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The magnetic contribution to the thermal conductivity /t mag , obtained from (1A8|) with 
-^ 22 ag nC3H . and the magnetic contribution to the resistivity p mag , following from (IA9D with 
P™ 3,6 (1C31) fulfill the Wiedemann-Franz law. It is, like in the case of scattering on the 
nonmagnetic impurities, the consequence of the relaxation time depending solely on the 
electron energy. 
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